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The set of Gaussian Type Orbitals g(n 1 ,n2,ng) of order (n + l)(n + 2)/2 and of common n = 
ni + 712 + «3 < 7, common center and exponential, is customized to define a set of 2n + 1 linear 
combinations t n ,m (~ n < m < n) such that each t n>m depends on the azimuthal and polar angle of 
the spherical coordinate system like the real or imaginary part of the associated Spherical Harmonic. 
(Results cover both Hermite and Cartesian Gaussian Type Orbitals.) Overlap, kinetic energy and 
Coulomb energy matrix elements are presented for generalized basis functions of the type r s t n , m 
(s — 0, 2, 4 . . .). In addition, normalization integrals J \g(ni, fi2, nz)\d 3 r are calculated up to n — 7 
and normalization integrals J \r s t n , m \d 3 r up to n = 5. 
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I. SCOPE 



z 2 . It is normalized as 56 



Gaussian Type Orbitals (GTO's) are widely used con- 
struction elements of basis sets in quantum chemistry. 
They became highly successful owing to fairly simple an- 
alytical representations of key integralsi»2iMi& that pay 
off in terms of speed when the matrix elements are calcu- 
lated. This work deals with the conversion of Cartesian 
or Hermite GTO's which have a product representation 
in Cartesian coordinates into harmonic oscillator func- 
tions which have a product representation in spherical 
coordinates. The latter are orthogonal with respect to 
two quantum indices — the total number of integrals to 
be calculated is reduced by the overlap and kinetic energy 
integrals over products of orbitals with the same center. 

Sec. |TT] introduces part of the notation. Sec. IIIII lists 
those linear combinations of Gaussians defined in Carte- 
sian coordinates that obtain the orthogonal property and 
are found to be the harmonic oscillator eigenfunctions. 
To extend coverage of the functional space, the oscillator 
functions are generalized in Sec.|lVl Section IVl computes 
two-center overlap, kinetic energy and Coulomb integrals 
of those without recourse to the expansions in Cartesian 
coordinates. Sec. IVII tabulates some of their products in 
support of 2-particle Coulomb Integrals. Absolute norms 
of orbitals play some role in Fitting Function techniques 
and are given in Sec. IVIII for the types of Gaussians dis- 
cussed before. 



II. TERMINOLOGY 



Let 



g(ni, 112,113) = a n/2 H ni (^/ax)H n2 (^y)H n3 (\/az)e ar 

(1) 

be a primitive Hermite GTO (HGTO) with exponent a 
and quantum number n = n\ + ri2 + 713 centered at the 
origin. H ni are Hermite Polynomials, and r 2 = x 2 + y 2 + 



3 

J 9 2 (n 1 ,n 2 ,n 3 )d 3 r=(^) 3/2 a n p(2n 3 -l)\\. (2) 

(ri + l)(n + 2)/2 different 3(^1,712)^3) exist for a given 
n. If n > 2, they build an overcomplete set of states 
compared to only 2n + 1 eigenstates Y,™ of the angular 
momentum operator. 

(— I < m < I) shall denote Spherical Harmonics in spher- 
ical coordinates^ and 



(4) 



generalized Legendre Polynomials. Their real- valued 
counterparts are 

y rn,c _ ^Y l m + Y l - m )/V2<xPr(cos6)cos(mip), (5) 
= -i(y™-y-" l )/V2cx P/"(cos0)sin(m^){6) 



and 



lf' c = If oc if (cos 0), 



(7) 



(0 < m < I), all normalized to unity, with parity (— ) , 
and orthogonal, 

smOdO f d^Yr> c (6, ^ m '' c (#, <p) 

sin9d6 f ^ dcpY l m ' s (e,cp)Y l f' s {d,ip) 
'0 Jo 

2tt 



(8) 



ammo I dipY™> c {B,y)Y™ ' s (6>,<p) = 0. (9) 
^0 

The main result of this work is support to use of GTO's 
in systems with heavy or highly polarized atoms by re- 
duction of the overcomplete sets to sets of 2n + 1 linearly 
independent combinations of GTO's. 
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III. COMPLETE SETS OF GAUSSIAN TYPE 
ORBITALS 

A. Hermite Basis 

Below, HGTO's are linearly combined into sets of real- 
valued functions t n>m (r,6,ip,a) (— n < m < n) that dis- 
play angular dependencies 4„, m oc Y™ c (9, ip) for m > 
and 4„, m oc Y~ m > s (6, ip) for m < by way of construction 
(see Appendix [A"f . Hence they are orthogonal 9 



r-2-K 

sinOdO / d(pt n , m t n ^ m i oc S n ^5 n 
o Jo 



(10) 



and complete with respect to the angular variables. Their 
norms arc listed in form of 



tl m d\ = N nm V2ir 3 a-*a n . 



(11) 



The N nm follow each time after three dots and are eas- 
ily derived from the expansion coefficients given and the 
overlap integrals 610 



g(n u n 2 ,n 3 )g(n' 1 ,n' 2 ,n' 3 )d r 



« ( " + "' )/2 (a 3/2 nH L ^ J+L ^ J K+-;-i) !! 



\2aJ 



(12) 

, all n j + rij even 
, any nj + n'j odd 



usingEq. (7.374.2) of Ref. OH or Eqs. (20) and (21) of 
Ref. LL2|. The radial dependence 4„. m oc (ar) n exp(— or 2 ) 
is found by inspection of the integrals mentioned in item 
1 of App. [XI which leads to 



Ln..m, — * 



(2n + 1) 



-{arTe- m - 2 Y™^ s \ (13) 



where 



Yl 



l{c . s} f Y™'°, rn>0 
Y~ m ' s , m<0 



(14) 



io, = 5(0,0 ) 0)...l/4; 
i li o=5(0,0,l)...l/4; 
t ltl = 5 (1, 0,0). ..1/4; 
t a ,o = 2g(0, 0, 2) - g(2, 0, 0) - g(0, 2, 0) ... 3; 
* 2 ,i =3(1,0,1)... 1/4; 
t 2 , 2 = 5(2,0,0) - 5 (0, 2,0)... 1; 
t2,_a=ff(l,l,0)...l/4; 
< 3 ,o = 23(0, 0, 3) - 3(3(0, 2, 1) + 5 (2, 0, 1)) . . . 15; 
t 8 ,i = 4 5 (1, 0, 2) - 3(3, 0, 0) - 3(1, 2,0)... 10; 
t 3l2 =3(2,0,l)-3(0,2,l)...l; 
i 3 ,- 2 = . 9 (1, 1,1)... 1/4; 



i 3 , 3 =3(3,0,0)-33(l,2,0)...6; 
U l0 - 83(0, 0, 4) + 63(2, 2, 0) - 24(3(2, 0, 2) + , 9 (0, 2, 2)) 
+3(3(4,0,0) +3(0, 4,0))... 1680; 
i 4 ,i = 43(1, 0, 3) - 33(1, 2, 1) - 3 3 (3, 0, 1) ... 42; 
£4,2 = 6(3(2, 0,2) -3(0, 2, 2)) 
-3(4, 0,0) + 3(0, 4,0)... 84; 

44, -2 = 63(1, 1, 2) - 3(1, 3, 0) - 3(3, 1,0)... 21; 

4 4 ,3=.9(3, 0,1) -33(1, 2,1)... 6; 
4 4 , 4 = ff(4, 0, 0) + 3(0, 4, 0) - 63(2, 2,0)... 48; 

4 4 , -4=.9(3,l,0) -3(1, 3,0)... 3; 

4 5 ,o = 83(0, 0, 5) + 15(3(4, 0, 1) + .9(0, 4, 1)) 
-40(3(2, 0, 3) + 3(0, 2, 3)) + 303(2, 2,1)... 15120; 
4 5>1 - .9(5, 0, 0) + 23(3, 2, 0) - 123(3, 0, 2) + 3(1,4, 0) 
+83(1,0,4)- 123(1, 2, 2)... 1008; 
4 5 ,2 = 2(3(2,0, 3) -3(0, 2, 3)) 
-3(4, 0,1) + 3(0, 4,1)... 36; 

4 5 , _ 2 = 23(1, 1, 3) - 3(3, 1, 1) - 3(1, 3, 1) ... 9; 
4 5 , 3 = 83(3, 0, 2) - 3(5, 0, 0) + 23(3, 2, 0) 

+33(1, 4,0) - 243(1, 2, 2)... 864; 
4 M = 3(4, 0, 1) + 3(0, 4, 1) - 63(2, 2,1)... 48; 

45, -4=3(3,1, 1)-3(1, 3,1). ..3; 

4 5 , 5 = 3(5, 0, 0) - 103(3, 2, 0) + 5.9(1, 4, 0) . . . 480; 
4 6 , = 163(0, 0, 6) - 5(3(0, 6, 0) + .9(6, 0, 0)) 
-120(3(0, 2, 4) + .9(2, 0, 4)) + 90(3(0, 4, 2) + 3(4, 0, 2)) 
+180 5 (2, 2, 2) - 15(3(2, 4, 0) + .9(4, 2, 0)) . . . 665280; 
4 6>1 = 83(1, 0, 5) - 203(1, 2, 3) + 53(1, 4, 1) 
-203(3, 0, 3) + 103(3, 2, 1) + 53(5, 0, 1) . . . 7920; 
4 6j2 = 3(6, 0,0) -3(0, 6,0) +3(4, 2,0) 
-3(2, 4,0) + 16(3(0, 4, 2) -3(4, 0,2)) 
+16(3(2, 0,4) -g(0, 2, 4))... 12672; 

4 6 , _ 2 = 163(1, 1, 4) - 16(3(1, 3, 2) + 3(3, 1, 2)) 
+3(1, 5, 0) + 3(5, 1,0) + 2.9(3, 3,0)... 3168; 
4 M = 83(3, 0, 3) - 243(1, 2, 3) + 93(1, 4, 1) 

+6.9(3, 2,1) -3.9(5, 0,1)... 3168; 
4 6 ,4 = 10(3(0, 4, 2) + 3(4, 0,2)) 
+5(3(2, 4,0) + 3(4, 2,0)) - 603(2, 2, 2) 
-3(0, 6,0) -3(6, 0,0)... 10560; 

4 6 ,_ 4 = 3(1,5,0)-3(5,1,0) 
+ 10(3(3, 1,2) -.9(1, 3, 2))... 660; 
4 6 , 5 = 3(5, 0, 1) + 53(1, 4, 1) - 103(3, 2,1)... 480; 
4 6 , 6 = 3(6,0,0) -3(0, 6,0) 
+ 15(3(2, 4,0) -3(4, 2,0))... 5760; 
4 6 ,_ 6 = 3(3(1, 5, 0) + 3(5, 1, 0)) - 103(3, 3,0)... 1440; 

4 7 , = 163(0, 0, 7) - 35(3(0, 6, 1) + 3(6, 0, 1)) 
-105(3(2, 4, 1) + 3(4, 2, 1)) - 168(3(0, 2, 5) + 3(2, 0, 5)) 
+210(3(0, 4, 3) + 3(4, 0, 3)) + 420 5 (2, 2,3)... 8648640; 
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t r ,i = 240 5 (3, 2, 2) - 5g(l, 6, 0) + 64 5 (1, 0, 6) 
-240.9(1, 2, 4) + 120g(l, 4, 2) - 240 5 (3, 0, 4) - 5 5 (7, 0, 0) 
-155(3, 4, 0) + 120^(5, 0, 2) - 155(5, 2,0)... 4942080; 

t 7 ,a = 15(5(6, 0,1) -5(0, 6,1)) 
+ 15(5(4, 2, 1) - 5(2, 4, 1)) + 48(5(2, 0, 5) - 5(0, 2, 5)) 
+80(5(0, 4, 3) - 5(4, 0, 3)) . . . 823680; 
i 7 ,_ 2 = 485(1, 1, 5) - 80(5(1, 3, 3) + 5(3, 1, 3)) 
+15(5(1, 5, 1) + 5(5, 1, 1)) + 305(3, 3,1)... 205920; 
t 7 ,3 - 35(7, 0, 0) - 2405(1, 2, 4) + 1805(1, 4, 2) 
-95(1, 6, 0) + 805(3, 0, 4) + 1205(3, 2, 2) 
-155(3, 4, 0) - 60 ff (5, 0, 2) - 35(5, 2,0)... 1647360; 
t 7A = 10(5(0, 4, 3) + , 9 (4, 0, 3)) - 3(5(6, 0, 1) + 5(0, 6, 1)) 
-60.g(2, 2, 3) + 15(5(2, 4, 1) - .9(4, 2, 1)) . . . 37440; 
t 7 ,_ 4 = 3(5(1, 5,1) -5(5, 1,1)) 
+ 10(5(3, 1,3) -.9(1, 3, 3))... 2340; 
t7, 8 = 60.9(1, 4, 2) - 5(7, 0, 0) - 55(1, 6, 0) - 120 3 (3, 2, 2) 
+5.9(3, 4, 0) + 125(5, 0, 2) + 95(5, 2,0)... 149760; 
t 7l8 = 0(6,0,1) -fl(0,6 s l) 
+ 15(5(2, 4,1) -,g(4, 2,1))... 5760; 
t 7 ,_6 = 3(5(1, 5, 1) + 5(5, 1, 1)) - 105(3, 3,1)... 1440; 
i 7l7 = 5(7, 0, 0) - 75(1, 6, 0) + 355(3, 4, 0) 
-21.9(5, 2,0)... 80640; 
i 8:0 = 128 5 (0, 0, 8) + 35(5(0, 8, 0) + 5(8, 0, 0)) 
+3360(5(0, 4, 4) + 5(4, 0, 4)) - 1792( 5 (0, 2, 6) + 5(2, 0, 6)) 

-1120(5(0, 6, 2) + .9(6, 0, 2)) + 67205(2, 2, 4) 
-3360(5(2, 4, 2) + .9(4, 2, 2)) + 140(5(2, 6, 0) + 5(6, 2, 0)) 
+210(4, 4,0)... 8302694400; 
t 8>8 = .9(8, 0, 0) + 5(0, 8, 0) - 28(5(6, 2, 0) + 5(2, 6, 0)) 
+705(4, 4,0)... 1290240; 
t 8 ,- s = g(7, 1,0) -.9(1,7,0) 
+7(5(3, 5,0) -5(5, 3,0))... 20160. 

t n ,-m with odd m are not shown explicitly, but are in- 
corporated implicitly by an interchange of the first two 
arguments of every 5 on the right hand side of i n , m and 
multiplication by (— 1)L TO / 2 J, like for example 

t 4 ,-i = 45(0, 1, 3) - 35(2, 1, 1) - 35(0, 3, 1) . . . 42. (15) 

This follows from applying the mirror operation x <-> y 
to the equations shown, which is <p <-> tt/2 — if in polar 
coordinates, and induces cos(rrup) «-» (— l)L m / 2 J sin(mtp) 
and Y r [ n - C <-► (-l)LW2jym, s if m is odd 
Application of the Laplace operator yields 



2a(2ar 2 - 2n - 3)t„ 



(16) 



The ij^m with exponential a = Mu)/(2h) are eigenfunc- 
tions of the 3-dimensional isotropic harmonic oscillator 



(- ^tt V 2 + ~Mco 2 r 2 )*„, m = (n + -)hujt n . m . (17) 



B. Cartesian Basis 

The previous list contains also the expansions in terms 
of Cartesian GTO's (CGTO's) as outlined in App. If 
each 5(711, 7i2, na) in the list is replaced by the primitive 
CGTO of the same triple index, 



f(ni,n 2 ,n 3 ) = x ni y n2 z n3 e~ 



(18) 



and each t n ^ m by t nim , the newly defined t n ^ m are cx 
r n exp(—ar 2 )Y™ c (m > 0) and cx r n exp(—ar 2 )Y r ^ ,n - s 
(m < 0), for example 



f 4 , 



4/(1, 0,3) -3/(1, 2,1) -3/(3, 0,1) 

4 /27T 



(19) 



The normalization integrals are recovered from the Kro- 
necker product of the vectors of the expansion coefficients 
and the overlap integrals^ 



f(ni,n 2 , n 3 )f(n' 1 ,n' 2 , n' 3 )d 3 r 

3 



I (JL) 

( 4Q! )(n+n')/2 V2a/ 



3/2 



IJ^+n^-l)! 



j'=i 



, all rij + even 
, any n 3 + n' odd. 



(20) 



to become 



II m d 3 r = N nm V2^a- 3 /{4a) r 



(21) 



[The similarity between Eq. JT2J) and Eq. (J20J) lets the 
iV nm already given in Eq. (|11|) show up again. The addi- 
tional sign in Eq. (I12p is positive if applied to 5-terms of 
a single t„ im , and does not mix things up.] 

A synopsis of Eqs. ([TTj) , and (|2l|) demonstrates 



tn,m — ^n,m/(2tt) • 



(22) 



An explicit expansion of the vector harmonics multi- 
plied with the Gaussian in terms of CGTO's is derived 
in App. EP£ 



r e-" ! r fcs) = _Lj22±i(n- \m\)l(n+\m\)\ 



cr 1 +cr 2 <l(n-\m\)/2\ 



E 



cri!u2! 



2tt 



\m\ 



CT1+CT2 



(|m| + 0-1 +<r 2 )! 



(n - \m\ - 2(ai + a 2 ))! 



(_)U/2J 



x /(H-i+ 2o "i.i + 2o "2,n-|m|-2(£ri + o-a)). (23) 

(In case of 1^™'° the sum over j attains only even values 
of j, in case of Y™' s only odd values.) This triple sum 
is t ni ± m up to a normalization factor that depends on n 
and m, and therefore (up to a different factor) also t ni ± m 
when / is replaced by 5. 
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C. Rayleigh-Type Expansion 

The expansion of a plane wave in terms of the t T 
reads 

oo / - l \ n 



n V * 

n 

E 



m— — n 



(2n + l)\\N n 
similar to the Rayleigh expansion 



-r n m ^ s >(k)i„, m (24) 



oo I 



e *" = 4tt£ E i l jl(kr)Yr(W m (i), (25) 



i=0 m=-i 



where r denotes the angular variables of the vector r. The 
expansion coefficients are the quotients of the integrals 
(e r \t ntm ) /{t ntm \t n>m ) . If r points into the z-direction, 
we have 



(26) 



whence 



All subsequent considerations assume that s is a non- 
negative even integer, (i) This ensures that t n ,m,s is an- 
other linear combination of HGTO's, which is made ex- 
plicit in App. [DJ Therefore a standard, indirect path 
of computing integrals is already established: linkage to 
known approaches'^! 4 - by dissociation of all £ n ,m,s i n 
the integrands. A relevant application is given in Sec. 
Vll (ii) The restriction leads to truncation of some series 
expressions that follow — important to numerical eval- 
uation — whereas odd s would sometimes yield infinite 
series, (iii) s + n > —3/2 is needed to guarantee con- 
vergence of the normalization integral, and (iv) the set 
of s = 0, 2, 4, . . . suffices to let the t n ^ m ^ s span the entire 
vector space of the GTO's. 

The generalization of Eq. (Til"]) reads 



tl m J 3 r = N nm + 2n „ + , 1) ; ! V2^^a n , (29) 



'(2n + l)!!(4a 
and the generalization of Eq. §TE\i 



1/3 
VH n , m , s (a,r) = s(2n+s + l)-j -4a I s + n+ - 



+4a 2 r 2 tn lW , a (a,r). 



2, 
(30) 



e i kz = y2e-^ 2 /(4 Q ) L + 2i * t _ 1 (!L)\ >0 

a 3 \a / 

i fk\ 3 1 fk\ 4 i fk\ 5 \ 

"15 W t3 '° + ^oUJ * 4 ' + 378oUJ * B '°" J 

The closely related Cartesian case is written down by 
replacing t n fi with (2a) n t n> o in Eq. (f2"7| for all n. 



t n ,m,s is an eigenfunction of a spherical potential with 
repulsive core as detailed in App. |EJ 



B. Fourier (Momentum) Representation 

The Fourier integral of the orbital centered at the ori- 
gin is related to Eq. (|B2[) . The Rayleigh expansion ([23)1 is 
inserted, then the radial integral is solved and re-written 
with Eqs. (11.4.28) and (13.1.27) of Ref.QJ: 



IV. GENERALIZED OSCILLATOR FUNCTIONS 



(a,k) = / e ikr t n , ro , s (a:,r)cfV 



A. Definition - a J 2 J Nn ' m (j1A n T ( n+ 2 + ^ 

V (2n+l)!! W r(n+|)a(3+ s )/ 2 



The basis of the t n>m defined in the previous sections -k 2 /(ia) p( s 3 k 2 y m { c > s }fl/Ri ) 

is incomplete with respect to the radial coordinate. The '' 1 2'"^" 2' 4a n 

combination g{2, 0, 0) + gr(0, 2, 0) + 5 (0, 0, 2) of HGTO's g ^ • ' 

contains a component oc r 2 e~ ar Yq ,c , for example, which (— s/2) a f^ 2 \ 

cannot be represented in terms of £ n>m . The subsequent ^ a\(n + 3/2)^ \4a/ 



sections therefore introduce a more general type of func- 
tions centered at space points A, which is derived from 2 / N nm ( s /2)! 
Eq. (JT3J) by multiplication with r s : = 4?r \j (2n + 1)!! q,(3+s)/2 

i n , m ,>,r-A) = 2"+y^^ iy |r-Ar(a|r-A|)« s l 2 UJ " W ' 

xe -a(r-A) 2 r ™{^4( r ^A). (28) where Pochhammer's Symbol 14 has been used. 
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V. INTEGRALS OF OSCILLATOR FUNCTIONS 
A. Overlap Integral 

The overlap integral is a convolution in real space and 
calculated with the shift theorem of Fourier analysis by 
use of Eq. (J3TJ) with t Jh „ hS (a,k) and t„/, m /, s ' (J3, -k): 



tn,m,a( a > T - A )*n',m',s' (A r ~ B)d 3 r 



+(")' 



Z n n' 



— m to' 



, >f C (C); 



m < m < : 

/ ^m+m' 



-mm 
i/irm' 



to + m —to —to / 1 



Z n n' 
m — to' —to to' 



y 



|m— m |,c 



tn,m, s (a,r)t„', m ', s '(/3,r - C)d 3 r 
d 3 fc 



V2 

to = to' < : 

-mm' = L ( I n n' 

yinri ^2 I to + to' — m —to' 



(2tt) 3 
= 2 n+ "'(47r) 



,(a,k)t 



Z n n 

— TO TO 



irN ntm N n < 



x 7 

s/2 



(2n + l)!!(2n' + l)!!a s + 3 /3 s '+ 3 

3 s' 
2+2 

s'/2 



- 7 C 2 



+ (-)' 

m > \m |, m < 

^ ^m+m' 



, ) y°' c (C); 



(35) 

|m+m' | , 

(C); (36) 

|m+m'|,c/Q\ 

(37) 



-mm _ 
hflnn' 



y/2 \m+ \m'\ —to — |to'| 



(n+n'+3)/2 r ( „ + 3 + * \ p ; _ + _ 



E 



2 2 



I 



(-s'/2) c 



7 

^ cr!r(n + | + a) W ^ cr'!r(n' + § + a') \/3 

n+n' „ , , 

53 (-)("-"'- ; )/ 2 ( 7 C 2 ) i / 2 r(^ + n + " +< + a + a') 



y TO — 1 77l' | —TO I TO' I 

< to < |to'|, to' < : 



I— |m' | , 



-mm _ 
c/lnn' 



V2 



I 

to + I to' 



n n 
-to — I to' 



(C); (38) 



Y m +\ m '\i s ^G') 



(-) 



Z=|n— n' | 

(n+?i'-/)/2+cr+o-' w n _ n ' 

x iCS' (6) g Hlxr 



z 



n n 

I TO— |m'| —TO I TO' 



y 



— m+|m | ,s 



(C). (39) 



-( 7 C : 



j r [Occurrences of Y^' s have to be interpreted as in Eq. 
(|3"9"| .] Some cases like to < 0, to' > or to' < to < are 
not covered by Eqs. (f33]l if39|) . but may be dispatched by 
( 32 )help of the permutation symmetries 



Here C = B — A and 7 = a/3/ (a + (3). The angu- 
lar integral y ; ™™ is expressed in terms of Wigner's 3j- 
Symbol^ (§3 3 >il (Chap. XIV) 



mm 
■nan 



)'{C)= yr(C) J sin9 k d9d<p k 

xYf(k)Y™^Hk)Y™ ' {c ' s} (k) 
(2Z + l)(2n+l)(2n' + l) ( I n n' \ _W/^ 

Eq. (4.6.3) of Ref. [l5j and adaptation to its sign conven- 
tion of Spherical Harmonics determine (C): 

to' = : 



Cn'(C)=C™(C) ; C?(C)=«(C). (40) 

In case of common center (C = 0), the orthogonality 
relation^ (§(4b)) 



i n , m , s (a,r- A)t n / m /, a / r - A)d 3 r 



-mm 
i/tan' 



(")' 

to > to' > : 



Z n n' 

TO —TO 



y 



m{c,s} 



(C); 



_ 2 2 "+ 3 ^ m ra r(n+| + f + ^) 

_VVro ' (2n+l)!! 7 (a + /3)(*+«'+3)/2 1 j 

results, which turns into Eq. (|29|) if a = /? and s = s' . [If 
the derivation starts from Eq. |32|) . the term J2a becomes 
l/r(/+3/2). The factor ( 7 c 2 )'/ 2 then means that nonzero 
overlap can occur only if the term I = is present in the 
sum over Z. In turn, I — becomes the only term that 
(33)contributes, and inspection of the lower limit Z = \n — n'\ 
introduces a factor S nn i . Then ?/orm = $mm' /(4tt), and 
eventually 



-mm 
yinn' 



V2 



(-y 



V2 

to = m > 



y 

to + to' —to —to' y ' 

Z n n' 
to — to' —to to' 



m+m ,c 



(C) 



(-«/2)«r /7 



s'/2 



y 



m — m ,c 



(C); 



(34) 



-fir v - 

a) \aJ 4^ cr' 



(-«72).' 



.'// 



y^2 y to + to' —to —to' 



y 



m+m ,c 



(C) 



s /2 

Jo ^!r(n + § + cr) \aJ ^ a'W(n + | + a') \0 
3 

xT(n + - + a + a') 

1 3 s s' 3 377. 



r(n + 4) 



2' 2 



2' a' /T 



6 



'/ 2 py'A> ri„ + 4) 



x 7 



n+5/2p 



r(n+f + f)r(n+§ + %) 



(42) 



2 2 



r ?i 



3 s' 



is calculated with Eqs. (9.180.2), (9.182.3) and (9.122.1) 
of Ref. [H] 



B. Kinetic Energy Integral 



s/2 

E 



(~*/2)a /7 



s'/2 



^a!r(n+ § + cr) W ^ a'!r(n + § + a') \/3 



E „/i 



(-*'/2)c 



X r(n+-+cr + cr' 



(45) 



Two "direct" methods exist to calculate the kinetic 
energy matrix. (i) Application of the kinetic en- 
ergy operator on t n > >m >, a >(P,r - B) yields three terms Coulomb integral 
oc t n ',-m',s'-2(/3,r - B), oc tn', m ',s'(fi,T - B) and oc 
t n ',m', s'+2{fl, r — B) by Eq. (|3"D|) if s' 7^ 0, and two terms 
by Eq. (fl6)) if s' = 0. Kinetic energy integrals thus reduce 
to two- or threefold application of Eq. (f3"2")) i 1918 common 
factors like y^™, could be reused, (ii) The second di- 
rect method is based on the conversion of the Laplace 
operator to a squared wave number in Fourier space: 



C. 1-Particle Coulomb Integral 

The Fourier technique also suits the attack on the 



t n , m , s (a,r- A)- 



1 



r — r 



-t, 



,((3,r-B)d 3 rd 3 r' 



= 4-7T 



d 3 k 



ikC 



(2tt) 3 k 2 

2 n+n' 167r 3 / 



~tn,m,8 (j^i k)tjy ,m' ,s' (/^; k) 



*n,m, a (a, r - A)V^„/, m /, s / 09, r - B)d 3 r 
d 3 k 



e- lk -^„, m , s (a, k)t n ,, m , )S ,(/3, -kp3) 



(2n + l)!!(2n' + l)!!a s + 3 f3 s '+ 3 ' 

^ + "' +1,/2r (»+H) r (»'+!+£) 



s/2 



("*/2)c 



s'/2 



The common result of both approaches may be summa- 
rized as 



E a!r(n+ | +a) 



n+n 



t„, m , s (a,r- A) ( — ^ ) i„/ im /, s /(/3,r-B)d 3 r 



-\n-\-n +1 



(4tt) 



x £ (_)(»-*'-0/2( 7 cr 2 )V 

/ — |n— n'| 

J — n — n. 



a) \a) cr'!r(n' 

a r(| + 2±2^ti+tr + (j / -l) 



r(z + §) 



n.m 1 y n' jrr 



_ 7 c- x C™(C)iFi(- 



(27i + l)!!(2n' + l)!!a s+3 /3 s '+ 3 



-<t-<t' + 1;(+-; 7 C 2 ). (46) 



The limit of common center is 



1 



s/2 

E 



(-s/2). /7 



s'/2 

E 3| 



(-s'/2)c 



7 

^ a!r(n + § + a) \a) ^ a>W(n> + § + a') \P 

n+n 

x 2 (_)(»-n'-0/2( 7C 2y/2 

£=|ro— ri' I 

_ /3 n + n' + l 

x r — 1 

1 2 2 



t n , mia (a,T- A)- 7T*n',m',s'(/?,r- A)d A rd 3 r' 

|r — r'| 

o2n+3 2 jy 
^ ly n,m 



(2n+ l)!!a( s + 3 )/ 2 /3( s '+ 3 )/ 2 



x 7 



3 s 



2 2 



3 s' 



a + a' + lj yS (C 



s/2 

E 



("*/2)c 



s'/2 



2 T - 



(s'/2) c 



o-!r(n + § + a) \aJ a'T(n + § + a') \/3 



(n+n'-l)/2+aW+l jl-n-7 
x V 1 ~ 



(7C 



2\CT 



(44) 



X r(n+-+cr + cr' 



(47) 



In the general case of I < n + n' or a + a 7 > 0, the 
which is Eq. (J32J) multiplied by 27 plus the replacement representation 



a + <r' — > a + tr' + 1 at three places in each of the last 
two lines of both equations. The limit of common center 
(C = 0) is 

t n ,m, s (a, r — A) f-^f) in',m',s'(/3,r- A)d 3 r 



a-a' + l;; + -; 7 C 2 ) 



(n+n'-i)/2+ff+</-l f i_ n _. 

E 



cr-^ + l)c 



£!(* + §) 



( 7 C 2 )-(48) 



2 2 "+ 4 7T7V n , r 



(2n + l)!!a( s+3 )/ 2 /3( s '+ 3 )/2 



lets Eq. gnj) look like Eq. ([32]) multiplied by 7r/7 plus 
the replacement a + a' — > cr + c' — 1 at three places. 
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The special case of I = n + n' = a + a' = may be 
paraphrased with Eq. (7.1.21) of Ref.fll 



1 F 1 (1;-, 7 C 2 ) 



7 C 2 2 



rfVTC 2 . (49) 



The remaining cases of I = n + n' > 0, er + a' = could 
be recursively attached to the representation (|49|) with 
Eq. (13.4.6) of Ref. H3, 

1^(1;* + ^7C 2 ) = {^(1^ + \, jc 2 ) i 

Some basic observations on the numerical evaluation 
of Eqs. dill), dH and (46]) are: 



Hi,. 



= if I + n + n' is odd. Therefore each second 



term in the sum over I may be skipped. 

• All appearances of the T-function have the form 
T(j + 3/2) = (2j + 1)!!V¥/2J' +1 with small integers 
j, and may be placed in a lookup table. (The asso- 
ciated factors \pK may be dropped, since the same 
number of T's appears in numerators and denomi- 
nators.) 

• In general, the number of terms in the sums is con- 
siderably smaller than the number of GTO pairs 
in the "dissociated" products. The focus shifts to 
efficient implementation of the yj%™, in Eqs. 



VI. TABLE OF PRODUCTS 



The 2-Particle Coulomb Integral 



|r — r' 



x *s,™, s -(/V - B)t n > tml<s >((3,T- B)tfVdV(51) 

may be computed by (i) expansion of the product 
tn,m,s{ot, r — A)t n i >m ',s'(f3, r — B) into a sum of terms oc 
t.,.,.(a+f3, r-P), all centered at P = (aA+j3B)/(a + /3), 
(ii) expansion of the product t Si ^|(a, r'— A)tn : m,s(P, r'— 
B) likewise into a sum of terms oc *.,.,. (a + /3,r — P), all 
centered at P = (a A. + (3B)/(a + /?), and (iii) comput- 
ing the remaining sum over 1-Particle Coulomb Integrals 
with centers at P and P as proposed in Sec. IV CI 

To implement the first two steps, an expansion of a 
product of two oscillator functions must be at hand. Its 
evaluation is straight forward as proposed here, though 
tedious: * n ,m,s(a ; r — A) and t n /. m ^ s >([3, r — B) are both 
expanded as described in App. [D] to obtain two sums 
of HGTO's. They are merged with the product rule of 
HGTO's: Those g are shifted to the common center P 



r — A = r — P 



B 



a + P 
a 

( 

a + (3 



(52) 
(53) 



where C = B A = (C x ,C y ,C z ). Combined at P, 
the product is a sum of g of quantum numbers up to 
n + n' + s + s', and finally re-converted into a sum over 
t.,.,.(o(+/3) r — P) as demonstrated in App.lAl To simplify 
the intermediate steps, one actually uses CGTO's and the 
associated 



*n 



T *n,m — *n,m,s/(2ck) 



(54) 



to work this out. Switching from the t representation to 
the / representation is done with the coefficients of the 
table in Sec. IIII Al To switch back one builds a table 
with the projection technique of App. |A] that starts as 
follows: 

/(0,G,0) = * ,o,o! 
/(o,o,i) = ti, ,o; 
/(0,l,0) = ti,_i, ; 
/(l,0,0) = t Ml0 ; 
/(0,0,2) = 1(^,0,0 + *D,o,a); 

/(0, 2, 0) = — £*2,0,0 — 5^2,2,0 + 3*0,0,2! 
f(2, 0, 0) = — g?2,0,0 + 5^2,2,0 + 3*0,0,2! 

/(0, 1,1) = £2,-1,0; 
/(l,Q,l) = i 2 , 1)0 ; 
/(l,l,0) = t 2 ,_ 2 ,o; 
/(0, 0, 3) = 5*3,0,0 + 1*1,0,2! 

/(0, 3, 0) = — 2*3,-3,0 — 20*3,-1,0 + 5*1,-1,2; 
/(0,l,2) = i(t3,_l >0 +*l,-l,2); 
/(1,0,2) = |(*3,1,0 + *1,1,2); 
/(0, 2, 1) = — j^*3,0,0 — 2"^ 3 ' 2 >° + iK 1 ' ' 2 ' 
/(2, 0, 1) = — j^y*3,0,0 + 2"^3,2,0 + 5*1,0,2! 
2, 0) = — j*3,3,0 — ^y*3,l,0 + 5*1,1,2; 
/(2, 1, 0) = 2*3,-3,0 — 20*3,-1,0 + g*l,-l,2; 

/(i,i,i) = t 3 ,_2,o; 

/(0,0,4) = gg-f^o.O + 5*0,0,4 + 7*2,0,2! 
/(2, 2, 0) = 2M*4,0,0 + T5*0,0,4 — ^ih,Q.2 ~ 8*4,4,0- 

Examples of products obtained with a Maple 2 ^ imple- 
mentation follow. The ~ signals the notational conven- 
tions that the first t on the left hand side has the ar- 
guments (a, r — A), that the second t on the left hand 
side has the arguments (/?, r — B), that all t on the right 
hand side have the arguments (a + /3, r — P), and that the 
right hand side is to be multiplied by exp(— 7C 2 ). The 
full writing of the first line would be 

£0,0,0 (a, r - A)t ,o,o(/3, r — B) 

= exp(- 7 C 2 )i ,o,o(a + /?,r-P). (55) 

Reduced exponentials a r = a/(a+/3) and j3 r = (3/(a+/3) 
are also inserted to simplify the notation. 

*o, 0,0*0, 0,0 — *o,o,o! 
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*0, 0,2*0, 0,0 — /3^C 2 <o,0,0 + *0,0,2 
*0,0, 2*0,0,2 — a r r C ?o,0,0 + *0,0,4 

+ {a 2 r + (3 2 - -a. r (3 r )C 2 i 0fia 

+2a r (3 r (a r — (3 r )C (C x ii t ifi + Cyt^-i^ + G^ti^o) 
— 2(a r — f3 r )(C x ti } i : 2 + Cyti-i^ + C z t\fi2) 
2 

+-a r (3 r (C - 3C z )t 2) o,o 

— 8a r f3 r (C x C y t2 : -2,o + C x C z t2Sfi + C y C z t2,-i,o) 
—2a r (3 r (C 2 — C 2 )t2,2,o] 
*i, o,o*o,o,o — f3rC z to t Q t Q + tififil 

*1, 1,0*0, 0,0 — f3 r G'xto,Q : + *l,l,0j 
*1, -1,0*0, 0,0 — PrCytofifi + tl,-l,o; 

2 - 2 ~ 

*1, -1,0*0,0,2 — a r l3 r C y C *0,0,0 + (A- — g«r)Cy*0,0,2 

+a r (a r C 2 - 2{3 r C y )ti-\fi + £1,-1,2 
—2a r f3 r (C x C y tii^o + Cj/Cz*i,o,o) 
— 2a r (C 2; i2,-2,o + C^*2,-i,o) 

+a r Gj,( -£2,0,0 + *2,2,o); 

*1, 0,0*0,0,2 — OL r /3 r C z C to, 0,0 + (A- — -^0! r )C z tofi,2 

+a r (a r C 2 - 2/3 r G 2 )f"i,o,o + *i,o,2 
—2a r (3 r (C x C z ti : ifi + C y C z ti t -i.o) 

2 

—2a r (C x t2,i,o + C y t2,-i,o) — gCt r C 2 t2,o,o; 



*1,-1,0*1,-1,0 



2~ l ~ 
-CtrPrC tofi,0 + 77*0,0,2 



— (a r — ^r)Cy*l,-l,0 — 77*2,0,0 — 7^*2, 2, 0! 
*1, 0,0*1, -1,0 — —0>rf3 r CyC z to,0,0 + *2,-l,0 

— a r Cy*i,o,o + PrGz*i,-i,o; 

*1, 0,0*1, 0,0 — ~ a rflrC 2 tofi_o + -(?2,0,0 + *0,0,2) 

— (a r - Pr)C z ti,o,o; 

*1, 1,0*1, -1,0 — —<%rPrC x C y tofi,Q + *2,-2,0 

— a r Cyti t i t o + PrGc*i,-i,o; 

*1, 1,0*1, 1,0 — —<%rPrC x to t Ofi + 7j*0,0,2 

— (a r — (3 r )C x ti 1 — 7*2,0,0 + 77*2,2,0- 

D 2 



VII. ABSOLUTE NORMS 

If GTO's represent orbitals or wave functions, their 
squares represent particle densities, and normalization 
follows from integrals like @ or (fTTj) . The first power 
rather than the second one specifies local densities, 
if these functions are the constituents of density fit- 



ting functions i 21 ' 22 i 23 The measure J g(n 1: n2,n 3 )d 3 r = 
(Tr/a) 3 / 2 5 n fi is generally zero and does not provide a use- 
ful substitute to Eqs. (fTTj) and ([2"9")) for that reason. 
The absolute norm as calculated below is the next sim- 
ple alternative. It quantifies how many particles have to 
be moved from some region of space to others to real- 
ize the specific density relocation, and becomes useful to 
qualify the relative importance of terms with products of 
expansion coefficients and the fitting functions. 

The integrated absolute value of 5(711, 7i2, TI3) is a prod- 
uct of three integrals, 



/ 



\g(ni,n 2 ,n 3 )\d 3 r = G(ni)G(n 2 )G(n 3 ) 



(56) 



with 



/oo 
a n > /2 \H n] (V^u)\e- au2 du 
-OO 

= 2oS n ^l 2 \ \H nj (u)\e~ u2 du. (57) 



G(rij) is a sum of 1 + [rij/2\ integrals delimited by the 
positive roots of H n . . Each integral is solved via Eq. 
(7.373.1) of Ref. [ill or Eq. (22.13.15) of Ref. OJ: 



G(0) 
G(l) 
G(2) 
G(3) 



= V^T" ~ 1. 772453850905516a' 1/2 ; 



= 2: 



4 v / 2a/e » 3.431055539842827a 1/2 ; 

4a (l + 4e~ 3/2 ) » 7.570082562374877a; 



G(4) = 4(a/e) 



3/2 



3 -V6/2 




G(5) 
G(6) 
G(7) 



3^/2^3 



19.855739152211958a 3/2 ; 
59.2575529009459587a 2 ; 
195.90006551027769a 5/2 ; 
704.821503307929489499a 3 . 



The integral J \f(ni, n 2 , n^)\d 3 r is calculated with ease 
via Eqs. (3.461.2) and (3.461.3) of Ref. HE 

The linear combinations defined in Sec. IIII Al may be 
evaluated in spherical coordinates and decompose into 
products of integrals over r [handled by Eq. (3.461) of 
Ref. [lE|, <p (yielding 2-7T or 4 for m — or m 7^ 0, re- 
spectively) and 9 [handled by determining the roots of 
P™(cos#) and decomposition into sub- intervals] . "Mo- 
nomic" cases like £1,1 , £2,1 or £3,-2, which relate t n>m to 
a single HGTO, are already represented by Eq. (|5o| and 
not listed again. 

/ 1*2 o\d 3 r = 8nJ— « 25.7190053432553290a' 1/2 ; 
J V 3a 
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*2,2 


d 3 r 


*3,0 


d 3 r 


*3,1 


d 3 r 


*3,2 


d 3 r 


*3,3 


d 3 r 


*4,0 


d 3 r 



1*4,1 
1*4,2 



1*4,- 





d 3 r 


j 1*4,4 


d 3 r 


f 1*4,-4 


d 3 r 


/ 1*5,0 


d 3 r 



8Vt7" ~ 14.1796308072441282a~ 1/2 ; read 



d 3 r 



d 3 r 



d 3 r 



104tt/5 « 65.345127194667699360; 

224 
~5~ 



77/1 

16arcsin(5~ 1/2 ) +4tt 

5 

49.94800887034627509; 
16; 

12tt w 37.69911184307751886; 

— 7rV357raf (V30 + 3) \/l5 - 2^30 
245 V v 



+(\/30- 3)^15 + 2^) 
765.99700145937577804a 1/2 ; 
8V^ (l + 128/7 3/2 ) 
112.180027342566280a 1/2 ; 
32V^ (l + 34/7 3/2 ) 

160.8439445477562992a 1/2 ; 
16^^ (l + 34/7 3/2 ) 

80.42197227387814961a 1/2 ; 

24V^ w 42.538892421732384655a 1/2 ; 

64 % /i^« 113.4370464579530257a 1/2 ; 

16V™ ~ 28.3592616144882564a 1/2 ; 



— 7T f 1701 + 640^70 
567 V 

2501.97074380428359a 



1*5,1 


d 3 r 


« 594.462027151417576a; 


1*5,2 


d 3 r 


= 1024a/9 = 113.7a; 


5,-2 


d 3 r 


= 512a/9 = 56.8a; 


1*5,3 


d 3 r 


= ^ (3064^2 - 2916 arcsin(l/3) 4 
81 V 






w 556.287108307150907a; 


*5,4 


d 3 r 


= 128a; 


5,-4 


d 3 r 


= 32a; 


*5,5 


d 3 r 


= 1207ra « 376.991118430775189a. 



y |*n, m , s (a, r - A)|d 3 r = a 



■a/2 



n + 3 



s/2 



|*n,m|^ ^* 

(58) 



VIII. SUMMARY 

A set of basis functions has been defined by recombina- 
tion of HGTO's or CGTO's such that the members with 
common center are orthogonal and complete with respect 
to the angular variables. They turn out to be cigcnfunc- 
tions of the isotropic harmonic oscillator. A specific gen- 
eralization of those allows (i) to keep the computational 
advantage of sparse overlap and kinetic energy integrals, 
(ii) backup by the tabulated GTO's in case of need for 
all integrals that are multilinear in the orbitals (2-particle 
Coulomb Integrals), still use of "direct" alternative for- 
mulas for overlap and kinetic energy integrals, and (iii) to 
maintain the vector space of functions covered by GTO's. 
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APPENDIX A: METHOD OF EXPANSION 



The absolute norms of the generalized oscillator functions 



The representations in Sec. IIII Al were obtained by re- 
engineering the straight-forward solution to the inverse 
problem as follows: 

1. expanding all (n + l)(n + 2)/2 GTO's of a fixed n 
into the complete set of Y™ (0 < I < n; —I < m < 
I). For each triple (rii, ria, 113), this expansion im- 
plies the calculation of (n+ 1) 2 integrals of the type 
/ g(n,i, ri2, n^Y™* sm. OdO dip . Most of these vanish 
due to the selection rules (i) ni+ri2+n 3 +l even and 
(ii) I + m + 713 even. One may also use that the ex- 
pansion coefficient in front of Y,~ m is the complex 
conjugate of the one in front of YJ m . The pedes- 
trian's way to calculate the integral over 9 and if 
is (a) transformation of g{n\ 1 712,^3) into spherical 
coordinates via Eq. (jTJ) and x — rcosc/?sin# etc. 
This results in a product of two integrals, one over 
9 and one over (p. (b) The integral over ip becomes 
elementary if all occurrences of sin ip and cos ip are 
substituted with Euler's formula, (c) The substi- 
tution cos 9 = u reduces the integral over 9 to a 
Hypergeometric Function 3F2 (p- 183 of Ref. [24]). 

2. gathering and recombining all Y" 1 within each ex- 
pansion in terms of y ; m,c and Y™' 8 . The list at 
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n < 3 (partial list at n = 4) then reads: 



and Yf" 1 ' 6 with 



-(-jLm/ajjrm,- ifmia even _ 



0,c. 



<?(0,0,0)~2Y C 
5(0,0,1)-^^ 
5(1,0,0)-^^ 



0, c 
L ! 

1, C 



5 (0,0,2)c(|^_4)y o ' c + -lf=i? 2 y 2 ' c ; 



'3" ^"0 1 3vV 

.9(1,0,1) c^K, 1 



?2 ^ c ; 



3(1,1,0)^^^ 

g (2, o, o) ~ (f i? 2 - 4)y °' c - ^ 2 y 2 °< c + ^i? 2 >f c ; 



0,c 



32 p3v°' c - 
5Vf 3 ' 



<?(0,0,3) ~ 1^(2^-5)^ 



g(l, 0, 2) ~ ^i?(2i? 2 - 5)1?'° + f sJl^Y^- 



5(1,1,1) 



16 ?jR 3y2, S . 



105 ' " 3 ' 
2 cwO,!: 16 T>3\r0,c 



g(2, 0, 1) ~ ^i?(2i? 2 - 5)Y X U ' C - ^R% 



16 
/105" 

2 



3v 2 > c - 



i? 3 K 



5 (2,l,0)c^i?(2i? 2 -5)r i 



rl.S 



i? 3 K 



3,s 



+8 



5(3,0,0) ~ § V3R(2R 2 - 5)Yl' c 



2 n3y3,c. 
35 2 3 ' 



£i? 3 K M 



5 (0, 0, 4) ~ (f R 4 - 32i? 2 + 24)y °' c 



^i? 2 (2i? 2 -7)y 2 u 



256 p4Y°' c - 
105 4 , 



ff (i,o, 3) ~ f J| (2i? 2 - 7)y 2 + iJ|i? 4 y 4 



rl,C 



(1,1,2)^^(2^-7)^ + ^ 
5 (2,0,2)^(f|i? 4 -fi? 2 + 8)y o f 



64 ?jR 4y2, s 



^0,c 



+ 27^ ( 2i? - 7 ) y 2 

128 e?4v°> c 



0,c 



16 



128 r>4vv 
105 2 4 



7V15 

64 n>4v 2 ' c 
2TV5 K ^ 



i? 2 (2i? 2 - 7)r 2 



2.c 



,(2, 1, 1) ~ ^i? 2 (2i? 2 - 7)Y 2 1 - S - §yflR 4 Y}> s 

rO,c , 32 n4 v 0,c 
r 1 a 



16 /Ip4y3,s. 
3 V 35 4 , 



32 R 2 
15"" 3 
32 D2/OR2 



5 (2,2,0)c(f|i? 4 -fi? 2 + 8)r u - 



105 J 



5(3,0,1)-- 



32 



3V35 U 4 ' 



|i? 2 (2i? 2 -7)r 2 1 



16 /2 n4ylc ,16 / 2 p4v 3 ' c 
-\/--tl J4 +TV 35 4 ' 



3. selecting the subset of equations that contain Y™' c 
for fixed m and n, yielding, say, 1 < q(m, n) < 
(n+ l)(n + 2)/2 equations. [One may generally find 
more stringent upper bounds for q by inspection of 
the selection rule (ii) given above.] The remaining 
task is to find q(m, n) numbers such that the linear 
combination of these q equations by these numbers 
does not contain any terms ^™„ c or Y["" ,s . This 
means computing a g-dimensional basis vector of 
a kernel of a matrix that contains all the expan- 
sion coefficients prior to those Y terms that are to 
be eliminated. Generally this matrix is non-square 
and r-dependent. 

4. normalizing this q-dimensional vector with some ar- 
bitrariness such that its components are small inte- 
gers and that t n ^ m /[(ar) n exp(— ar 2 )Y™> c ] are pos- 
itive numbers. Those q components are the expan- 
sion coefficients in front of the g(n±, 112,713) of the 
table for m > 0. Counting terms shows q(6, 2) = 8, 
q(6,0) = 10 and q(3,2) = 2, for example. Experi- 
ence shows that — up to this normalization factor 

— the expansions are unique for n < 8 at least, i.e., 
the aforementioned kernel is one-dimensional. 

5. performing steps 3 to 4 for the Y™' s in an equiv- 
alent manner by elimination of terms oc Y^^ and 
oc Y[ a ' c to obtain t nt - m . 

An additional shortcut exists once t n _ m and t ni _ m are 
known for a specific "anchorage" m. Application of the 
ladder operators 

L± = L x ± iL y = i(yd z - zd y ) ± (xd z - zd x ) (Al) 

of angular momentum quantum mechanics on Y™ yields 
l^ 71 . Decomposition of the two equations 

-^±(^71, m ~l~ ^n,— m) ^ ^n,m±l ~l~ ^n, — m=pl (-^-^) 

into 4 real- valued equations yields a similar recursion for 
t n , m . Effectively one applies 

(yd z - zd y )g(n 1 ,n 2 ,n 3 ) 

= n 3 g(n 1 ,n 2 + l,n 3 - 1) - n 2 g(ni,n 2 - l,n 3 + 1); 
(xd z - zd x )g(ni,n2,n 3 ) 

= n 3 g{ni + 1, n 2 , n 3 - 1) - nig{m - 1, n 2 , n 3 + 1) 

term by term to a pair of equally normalized t niTn and 

The quickest alternative to obtain the expansion is de- 
veloped in App. [Cl further down. 



The symbol ~ means a factor e~ ar a n l 2 ^pK has 
been omitted each time at the right hand side, 
and R stands short for -Jar. The coefficients of 
g(n 2 ,ni,n 3 ) are derived from those of g(nx, n%, 713) 
by replacing Y™' c with (— ) L m / 2 J y™' c if m is even, 
Y, m ' c with (— ) L m / 2 J y ; m ' s anc [ v j ce ve rsa if m is odd, 



APPENDIX B: CORRESPONDENCE BETWEEN 
CARTESIAN AND HERMITE GTO'S 

A Fourier Transform switches from Cartesian to Her- 
mite GTO's and vice versa^ 6 . All HGTO's g(n u n 2 , n 3 ) 
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at the right hand side of the equations of Sec. QTT3 trans- 
form into their associated CGTO's /(ni,n 2 ,n 3 ) in k- 
space, 

[ e ik -*g(ni,n2,n 3 )d 3 r=(l) 3/ % n e- k2 /wf[k?, 

3=1 

(Bl) 

and the components within the t n , m (r,O r ,tp r ,oi) at the 
left hand sides keep their angular dependence, 

e*- r (ar) n e~ ar2 Y™(e r ,^ r )d 3 r 

= © 3 ^- fc2/(4Q) V(^)- (B2) 

Subsequent re-insertion of r for k leads to the expansion 
of y„ m ' c or F„ m ' s in terms of CGTO's. Residual dangling 
factors depend on n, but not on individual rij, and are 
constant within each expansion, which allows to stay with 
the coefficients as shown. 

APPENDIX C: DIRECT EXPANSION FORMULA 

The observation in App. [B] leads to a closed formula of 
the linear expansions of Sec. IIII A\ because it suffices to 
expand t n<m and actually r n Y™^ c ' s ^ (6, <p)e~ ar in terms 
of /(ni, 7i2, 713). An explicit writing of Eq. H) is given 
by a re-formulation of Eq. (8.812) in Ref.[rp 

if (cos 0) = 2 l sm m 9 



Whence 



Vl'm+v + 1 — l\ cos" 6 



u=0 



Terms with odd I — m — v do not contribute to the sum 
(ICTI) . A synopsis of Eqs. (1.331.1-3) ofRef.GHis 



cos(mip) 
sin(mip) 



E (- 



1 LJ/2J 



• = r 0,3,4.. 

■J 1 1,3,5.. 



cos m 1 (ps'm J <p. 



(C2) 

The product of Eqs. (|C1[) and (|C2|) is inserted into Eqs. 
([5]) and ©. Simple counts of the powers of the trigono- 
metric functions constitute^ 



m + v + 1 — n 



n n— v— m 



v\{n — m — v)\ 



m 



x £ (-)ti/ 2 H~:) x m-3 y 3 zr (C3 ) 
J=0(1) VJ 



(The sum is over even j for Y™' c and over odd j for 
y^™' 5 .) The residual even power of r is also broken down 
to Cartesian components 



ri — v—ra 



E 



n — v—m \ 1 



(C4) 



(n-|m|)/2 



tn.m 0^ 



E 



i/=0 
Bx+02<V |m| 

x E E 



(2n-2i/- l)!!(-) w 
2 ly (n - |m| - 2z/)! 

(_)Lj/2J 



<7i!cr 2 !(z/ - cr x - cr 2 )! \ j 



0-1,0-2=0 j=0(l) 

x /(M - j + 2cti, j + 2ct 2 , n- |m| - 2(<ti + cr 2 )). 

(The sum over j is restricted to even j if m > 0, and odd 
j if m < 0.) Resummation enables performing the sum 
over v analytically to obtain Eq. (|23|) . 



APPENDIX D: EXPANSION OF GENERALIZED 
ORBITALS 

The expansion of t n>m , s — r s t n ^ m in terms of HGTO's 
is obtained as follows: i„, m is expanded as listed in Sec. 
IIII A( each g(n\,ni-,n%) within the sum is written as a 
product according to Eq. |T]). r s — (x 2 + y 2 + z 2 ) s l 2 is 
expanded in Cartesian coordinates as in Eq. (|C4[) . and 
each of its Cartesian components recombined with the 
Hermite Polynomial of the same Cartesian direction by 
repeated application of 

2xH ni (x) = H ni+1 (x) + 2nH ni -x{x). (Dl) 

s = 2 just needs 



x 2 H ni (x) 



\H ni+2 (x) + {^ + l)H ni {x) 

+n 1 (ni-l)H ni _ 2 (x), (ni > 2); 
±ff 3 (aO + §ffi(aO, 



\H 2 (x) + lH (x), 



(»i = 1); 
(m = 0). 

(D2) 

and the equivalent expressions for y 2 H ri2 (y) and 
z 2 H n3 (z). The powers s = 4, 6, . . . follow by recursion 
of Eq. (ID21 . 

The case of CGTO's, t n ,m,s = r s t n , mi proceeds even 
simpler, because the obvious x s f(n\ 1 n 2 , TI3) — f{ni + 
s, 77.2,71.3) etc. does not need auxiliary formulas like Eq. 
(ID2p to merge the components of r 2 with the / (ni, n 2) ^3) 
contained in Sec. IIII At 



APPENDIX E: OSCILLATOR WITH REPULSIVE 
CORE 

As Eq. (j3"0|) suggests, t n<niiS is an eigenfunction of the 
spherical harmonic oscillator with superimposed repul- 
sive core potential 



(El) 



The standard discussion of the single-particle 
Schrodinger equation deals with a fixed potential, 
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hence fixed Q. Its eigenf unctions ip nm j are found with 
the method proposed in Ref. H3: 

^■(a,r) = e — >-i 1 F 1 (-j;p+ \; 2ar 2 )Y^ {c - s} (r) , 

(E2) 



where 



p = - + J\ + Q + n{n + l), 



(E3) 



and j = 0, 1, 2, . . . is the number of nodes in the radial 
function (the zero at r = not counted) . The 2n + 1-fold 
degenerate energy eigenvalues are 



E 



1. 



J -—Aa {]+P+ -). 



(E4) 



[Additional "accidental" degeneracies occur if p is in- 
teger, which includes all cases of even, integer Q = 
(q + n)(q — n—l) up to 200 for example.] The Confluent 
Hypergeometric Functions may be written as Laguerre 
Polynomials 



1 



iFt(-tv+ --2ar 2 ) - - L { ?~ X,2) 



(E5) 



which helps to derive the orthogonality relation with Eq. 
(7.414.3) of Ref. IH 



2(2a)P+ 1 /2 - 
(E6) 

The t n ,m,s defined in Eq. (|28p refers to the case of integer 
Q = s(s + 2n + 1) , integer p = s + n + 1, and j '• = 0. Since 
the factor r s (ar) n in (f2"5| may be written as a linear 
combination of the polynomials r p ~ l iFi() in Eq. (|E2[) 
and vice versa, the set of basis functions t n , mtS (a, r) spans 

the same phase space as the set of ^nrnji a ^ Y )- 

A computational advantage of the eigenfunctions (|E2p 
is the orthogonality with respect to three indices man- 
ifested in Eq. (|E6|) . whereas the t n ^ m , s are orthogonal 
with respect to n and m but not with respect to s [see 
Eq. (|4"l"j) ]. — consequence of the fact that the if> n mj are 
eigenfunctions of a fixed potential and Hermite operator, 
whereas the various t n>m _ s belong to different potentials 
(different Q). A disadvantage of the eigenfunctions (|E2|) 
is that their Fourier representation has rather compli- 
cated sets of polynomials which replace i-FiQ in Eq. (|3ip . 
which takes negative effect on the complexity of all for- 
mulas derived in Sec. W\ 
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